Langmuir wave linear evolution in inhomogeneous nonstationary anisotropic plasma 
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Equations describing the linear evolution of a non-dissipative Langmuir wave in inhomogeneous 
nonstationary anisotropic plasma without magnetic field are derived in the geometrical optics ap- 
proximation. A continuity equation is obtained for the wave action density, and the conditions for 
the action conservation are formulated. In homogeneous plasma, the wave field E universally scales 
with the electron density N as E oc iV 3//4 , whereas the wavevector evolution varies depending on 
the wave geometry. 
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I. INTRODUCTION 

The energy of a wave propagating in nonstationary 
medium can be manipulated by controlling how the pa- 
rameters of the medium evolve 0, [11 IS, 0, IE H 0] > f° r 
example, the energy can be pumped up, transported, fo- 
cused, and (or) deposited where necessary. In the case 
of Langmuir, or plasma waves [![, there may be im- 
portant high energy density applications connected with 
compressing plasma targets, because these targets may 
advertently or inadvertently contain wave packets that 
would be amplified along with the densification. Hence 
understanding of the Langmuir wave evolution in nonsta- 
tionary plasma is needed. To develop such understanding 
is the purpose of this paper. 

Specifically, we assume the geometrical optics (GO) 
limit 0, E3, El, E2, EH, when the plasma parameters 
vary sufficiently slowly in time and space. We also as- 
sume that a wave is linear [l4| , and no collisions, ioniza- 
tion, or recombination take place [l5| . In this case, the 
plasma dynamics should allow a Lagrangian formulation 
[i~6l [Tm] ; thus it is anticipated to comply with the gen- 
eral theorem of GO which states that the wave action 
is conserved in inhomogeneous nonstationary medium 
[H, [H, [H SI IH HI H. Previously, the theorem 
was independently rederived for a variety of oscillations 
[3 El, HE HE H3, HI , confirming the general treatment; 
particularly, space-charge waves in cold electron beams 
were considered, similar to Langmuir waves in cold plas- 
mas [HHO, [M[. However, for thermal plasmas there has 
been less agreement, and some of the models proposed in 
literature do not comply with the action conservation. 

The Langmuir wave action, or "plasmon" conservation 
theorem (PCT) was reported in Refs. [H, [H, [H], HU , ac- 
counting for nonlinear effects; however, the inhomogene- 
ity of the background plasma was neglected there (see 
also Ref. (111). The density inhomogeneity was included 
in a linear treatment in Ref. [37j . yet within a model 
assuming constant temperature. More precise models 
of Langmuir waves in inhomogeneous plasma (see, e.g., 
Ref. [38| and references therein) did not specifically ad- 
dress PCT and assumed stationary medium; also, the 
wave equation derived in Ref. [3!| is not entirely correct 
(see, e.g., Refs. [io|, 5l[ and Sec. IIV Ap and hence is at 



variance with the theorem. Similarly, the kinetic mod- 
els offered in Refs. [H |43[ are erroneous, as explained 
in Refs. [3, 5f|, and so is the corresponding part of 
Ref. as argued in Appendix [Bj Thus, the accuracy 
of PCT with respect to the temperature corrections was 
not fully assessed. 



An accurate kinetic treatment of the thermal effects 
was eventually proposed in Refs. 0, [45j]. Particularly, 
it was shown that the Langmuir wave action in isotropic 
nonstationary collisionlcss plasma is conserved in the GO 
limit, assuming that Landau damping is insignificant. 
However, the solution in Refs. 0, [45[ is incomplete, be- 
cause collisionless plasma may not remain isotropic in 
the presence of inhomogeneous average flow [471 ]. Thus 
it yet remains to derive an explicit equation for a Lang- 
muir wave in nonstationary inhomogeneous anisotropic 
plasma and show how the wave parameters evolve. 



These results are reported in the present paper, which 
thus completes the studies in Refs. [H [3?], 53, EH and 
finally reconciles the Langmuir wave dynamics in inho- 
mogeneous nonstationary warm plasmas with the general 
principles of the Lagrangian GO. Specifically, we derive 
a continuity equation for the wave action density and the 
explicit conditions under which the action is conserved. 
Hence it is shown that, in homogeneous plasma carrying a 
Langmuir wave, the wave field universally scales with the 
electron density N as E oc iV 3 / 4 , whereas the wavevector 
evolution varies depending on the wave geometry. 



The paper is organized as follows. In Sec. [IT] we intro- 
duce our basic equations. In Sec. Mil we find the Lang- 
muir wave dispersion relation in homogeneous anisotropic 
plasma. In Sec. IIVI we derive the equations for GO rays 
and the amplitude of Langmuir oscillations in inhomo- 
geneous nonstationary plasma. In Sec. [V] we use those 
to obtain a continuity equation for the wave action den- 
sity; we also derive the scalings for the oscillation field 
amplitude and wavenumber. In Sec. I VII we summarize 
our main results. Supplementary calculations are given 
in appendixes. 
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II. BASIC EQUATIONS 



III. HOMOGENEOUS STATIONARY PLASMA 



Consider a Langmuir wave in unmagnetizcd nonrela- 
tivistic plasma with given flow velocity V(r, t). Neglect 
ion oscillations and assume collisions and Landau damp- 
ing to be insignificant on time scales of interest. For 
electrons adopt the low-temperature approximation re- 
quiring /cAd <C 1, where k is the wavenumber, and Ad is 
the Debye length (cf. Eq. (11) in Ref. HI). Hence an 
asymptotic closure of the hydrodynamic model is pos- 
sible, via omitting the heat flux, and one obtains, by 
taking the first three velocity moments of the electron 
Vlasov equation 0, [H, H : 



3 t N e + V • (iV e V e ) = 0, 



dtV e + (V e ■ V)V e 



E v^-— ^ 

m e N e m e 



(1) 
(2) 



d t P e + (V e ■ V)P e + P e (V • V e ) + 

+ [(P e V)V e ] + [(P e V)V e ] T = 0. (3) 

Here N e is the electron density, V e is the electron flow 
velocity, e < and m e are the electron charge and mass, 
P e is the electron pressure tensor (which is symmetric 
by definition), E is the average electric field (if any), ip 
is the wave electrostatic potential [slJ , and the index T 
denotes transposition. Separate the slow and the quiver 
variables, correspondingly, as 

N e = N + N, V e =V + V, P e =P + P, (4) 

and assume that the oscillations are weak, i.e., 

n = N/N < 1, (5) 

and similarly for the pressure. (However, the absolute 
value of V is unimportant.) Then the following linear 
equations are obtained: 

<9' f n + h-V + V-V = 0, (6) 

d' t V + (V • V)V - n + (e/m e ) = 0, (7) 

8'tP + (V • V)P + P (TrW) + (WP) + (WP) T + 

+ P(V-V) + [(PV)V] + [(PV)V] T = 0, (8) 

n = (nV • P - V • P)/(m e N), (9) 

VV = -K/e)^n. (10) 

Here we introduced the partial time derivative in the 
frame of reference K' (further denoted by prime) moving 
with velocity V with respect to the laboratory frame K: 



d' t = d t + (V- V). 



(11) 



In addition, we introduced h = V IniV, W = VV (which 
is a tensor with elements Wjt = dVj/dxg), the electron 
plasma frequency u> p — *J4ire 2 N/m e , and Tr for the 
tensor trace. The complex notation is also henceforth 
assumed for the quiver variables. 



In the case of homogeneous plasma with V = const, 
the exact eigenmodes of Eqs. (pj)- (Ti"U|) are found as fol- 
lows. Assume 

n, V, P ~exp(ik-r-iwi). (12) 
Then one gets 

-iuj'n + ik-V = 0, (13) 
-io/V + ik(e/m e ) <p - II = 0, (14) 
-iu'P + iP(k • V) + i[(Pk)V] + i[(Pk)V] T = 0, (15) 

n = -ikP/(m e N), (16) 
-k 2 ip — -(m e /e) u>pU, (17) 



with i 



lu — k-V being the wave frequency in the frame 



where the plasma average flow rests. Hence 
HV = cj' 2 V, 



(18) 



i.e., V must be an eigenvector of the tensor H given by 
1 



H 



kk 2 



2(Tk)k+ (k- Tk)l 



(19) 



where kk is a dyad, 1 is a unit tensor, and T = P/N is 
the electron unperturbed temperature tensor. 

For a longitudinal wave one has from Eq. (TT3]) that 



V = ncu'k/k 2 . 
Substituting this into Eq. fTS]) yields 

, 2 # , (k • f k) 



Tk 



2P 



(20) 



(21) 



Since the coefficient on the right-hand side is a scalar, k 
must be an eigenvector of T: 



Tk = m e w^ e k, 



(22) 



where the eigenvalues satisfy v\ e > 0, because T is pos- 
itive defined; thus v^e can then be understood as the 
electron thermal speed along k. Then 



ft(k)+k- V, 



(23) 



where the function fi(k) (numerically equal to lu') is 
given by fi 2 = uJi 



3fc 2 w 2 e (cf., e.g., Ref. [|5J, Chap. 8], 



Refs. [53j, [5J, [55| for isotropic plasma), or 



n 2 = 



k-Ck, C = 3T/m e 



(24) 



Therefore, similarly to elastic media [56j,l57J,l58, 59, 60], 
a longitudinal wave in plasma must have k along a prin- 
cipal axis of the temperature tensor [6l[ . Since T is sym- 
metric, there always exist at least three such directions. 
On the other hand, in isotropic plasma, any axis can be 
considered a principle axis of T; then Langmuir waves 
can propagate in arbitrary direction. 
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IV. GEOMETRICAL OPTICS EQUATIONS 
A. Approximate wave equation 

Consider now a more general case of inhomogeneous 
nonstationary plasma. Assume, however, that the wave 
amplitude, the frequency, and local wavevector vary on 
large temporal and spatial scales T and L in the plasma 
average flow rest frame K' . Specifically, we assume 



e = l/min{w'T,fc'L} < f, 



(25) 



with k' = k in the nonrelativistic limit used here. In 
this case, henceforth referred as the geometrical optics 
(GO) limit 0, [Tol | , an approximate scalar equation for a 
Langmuir wave can be obtained which will capture effects 
to the leading order in the parameter e. 

To derive this equation, first, apply d' t to Eq. © to get 



d' 2 n + d' t (h ■ V) + d' t (V • V) = 0. 



(26) 



The third term on the left-hand side here is found by 
taking the divergence of Eq. (J7]) : 

d't (V ■ V) = co 2 p n + v • n - (V • V)(V • V) - {V, V}, 
where we substituted Eq. (fT0|) and introduced 



{V, V} = V ■ [(V • V)V + (V ■ V)V] - 

- (V • V)(V • V) - (V • V)(V • V), (27) 

accounting for the fact that d't and V may not commute. 
Therefore Eq. (|26[) takes the form 

d' 2 n + Lo 2 p n + V • II + R = 0, (28) 



with R given by 

R = -(V ■ V)(V • V) + d't (h • V) - {V, V}. (29) 

The first term in Eq. ([29]) is of order e 2 and can be 
neglected. The second term is evaluated as 



d'(h ■ V) ^ d'n n k 

^^" 2 5F^ V ^'P (30) 



(d'/dt = d't), where we used that, to the zeroth order 
in e, one can employ Eq. (|20|) and 

d't w -ifi, V w ik. (31) 

The third term in Eq. (|2T)|) can be put as (Appendix [XJ 

{V, V} = 2(Wj) • (Wj) - 2(V x V) • (V x V), (32) 

where we assume summation over repeated indexes. Both 
curls are of order e here; thus {V, V} w 2(VVj) ■ (Wj). 
Eqs. 0, (EH) further yield 



{V, V} w -2 (a' t n) kW • k/fc 2 . 
Hence Eq. ([2^)1 rewrites as 



R = 2 



9'n / 



Vw„ + kW] • 



fc 2 ' 



(33) 



(34) 



An approximate expression for V • II is obtained simi- 
larly ( Appendix |B|) . Substituting that and Eq. (|3"4")l into 
Eq. ([28]) one gets for isotropic plasma [621 ] : 



d ,2 n 
9i 2 



w 2 n - 3v^ e S7 2 n + 2 



9'n / Q 
dt 



kW| 



and in the general case of anisotropic plasma, which we study below, 



d' 



dt 2 



+ w 2 n - Cji 



d 2 



d'n ( duj. 



dxj dxe 



k 2 



6Vn • Vu^e = °> 



k.jk„\ dC,f dn 



k 2 J dxj dxi 



(35) 



(36) 



B. Eikonal equation 

Equation (f3"6"]) can be solved using the GO approach 
d, HI], specifically as follows. Take 



n = Me l \ 



where Af is the slowly varying envelope. Substitute 
Eq. (|37p into Eq. (|36|) and first consider the terms of 



order e°; hence the eikonal equation 

[-{d t e - we ■ v) 2 + uj 2 v + ve ■ eve] jv = o. 

(37) Since, by definition, 

d t 9 = -lu, V6» = k, 



(38) 



(39) 
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Eq. (I38|) is equivalent to Eq. ([23]) , except now the plasma 
parameters may slowly depend on r and t: 



= n(k;r,t)+k-V(r,t). 



(40) 



Differentiate Eq. (|40p with respect to t and with respect 
to r and use Vw = — (?tk, flowing from Eqs. (I39)) . Then 

dtw = 9tw(k,r,t), d t k = -Vw(k, r, i), (41) 

where the partial derivatives are taken at fixed k; also 

dt = d t + (v s ■ V), (42) 

and v g = 9kw(k, r, i) is the group velocity: 



U + V, U = <9 k O(k,r,i). 



(43) 



Since v „ eq uals the velocity at which the envelope prop- 
agates [|52j, Chap. 4], one can also write 



d t r = <9kW. 



(44) 



Together, Eqs. (|4Tj) . (|44j) are known as GO ray equations 
[|52j, Chap. 4] and can be considered as canonical equa- 
tions [with the Hamiltonian hLo(k,r,t); Eq. ([40]) ] which 
determine the dynamics of "plasmons", i.e., quasiparti- 
cles with velocity v s , momentum Kk, and energy huj; see 
also Ref. [H HHf. 



C. Amplitude equation 



The equation obtained from Eq. 
in e reads 



in the first order 



where k is a function of r and t [unlike in Eqs. (|41[) , where k is an independent variable], and 

The same expression can be written also as follows. Use Eq. (|4Tj) for o? t k to get 

dfcf uj p duj p kjk s dC js 

-dt=-i[d^-xT^- kjWjt - (47) 

Hence Eq. ([4"u| is put in the form 

_ _ k^ dk^ / fi^ _ t^p\ fc^ _ fcjfcifcs , 

k 2 dt \u) p n J k 2 dx t 2ilk 2 dx e ' 1 ' 

We now use the expression for f2 [Eq. 1[2"4"[)]. Vw p /w p = VN/(2N), and 

k • dtk/fc = dtife. (49) 

Hence Eq. (|4l)]) can be represented as 

r = -^ + ^( C ^-^J- (50) 

Since the wave is assumed propagating along a local principal axis of the temperature tensor (64| . one also has 
Cjtkjki/k 2 = 3w^ e and 3v^ e k s = C s tki\ thus Eq. (|5U|) can be further put as 

rflnfc t/j-fcj k 3 k e k s dCj£ 
dt 2 ~ 20fc 2 9x s ' [ ' 

where we used Uj = Cj£k f /Q [Eq. (|4*3"|) ]. Then Eq. (|4"5|) rewrites as 

a'i(0|AT| 2 ) + V • (n\Af\ 2 U) + n\Af\ 2 (U • h - dt lnfc 2 ) = 0. (52) 
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We now use 



d' t (n\Af\ 2 ) + V • (ft|A/] 2 U) = d t (n\AT\ 2 ) + n\Af\ 2 \7 • U, U • h = d t hxN + V ■ V, 



r 



(53) 



the latter being due to 

d' t N = -NV ■ V, (54) 

which flows from Eq. {J). Hence Eq. ([52]) reads 

d t \n(Q\Af\ 2 ) + d t In N - d t In k 2 + V • v g = 0. (55) 

In principle, this allows one to calculate the envelope am- 
plitude | TV" | along the GO rays, as discussed in Sec. IV Al 

V. DISCUSSION 

A. Wave action. Equation of state 

Introduce the wave average energy density 
\E'\ 2 d(e[u>) 



£' 



Wtt dbj' 



(56) 



in the frame K' traveling with velocity V. Here E « E, 
and the longitudinal permittivity in K 1 equals that in 
the laboratory frame K: e\ = £\ [65j. Neglecting the 
corrections due to finite e and using Eqs. ©-(HU), £i is 
derived like for isotropic plasma [[521]. Chap. 3]: 



ei = l- 



3k 2 v 2 e 



(57) 



(Here k parallel to the principal axis of T is assumed, as 
before [64]].) Thus Eq. ([56]) rewrites as 



£' 



n 2 \e\ 2 



w 2 8ir 



(58) 



[where we used ei(fl + k • V; k) =0], so this energy den- 
sity is always positive, unlike that in K [66] ]. 

Further, define the wave action density J, or the num- 
ber of quanta (plasmons) per unit volume, as J = £' /uj' 
[H,|67|, where we take uj' > 0, by analogy with discrete 
systems (see, e.g., Sec. Ill of Ref. (68[). Then 



J oc n\AT\ 2 N/h 2 



(59) 



where we used E 



-i\fnp and Eq. (1 1 7[) for ip. From 
Eq. ([55]) , it follows then that d t In J + V • v g = 0, or 

d t J + JV-v 3 = 0. (60) 

The latter is also equivalent to a continuity equation: 



a f j + v-(v s j) = o. 



(61) 



Hence the wave total action is conserved, 

J d 3 r = inv, (62) 

and the dynamics is thereby called adiabatic. 

Eqs. ([60 ]) -([62 ]) agree with the previous results for 
space-charge waves in cold plasmas and electron beams 
[29l [30l [3l| , as well as phenomenological hydrodynamical 
treatment of the corrections due to the electron homoge- 
neous temperature (37| and kinetic treatment for inho- 
mogeneous nonstationary but isotropic plasmas [44l |45| . 
By construction [H, [13], the hydrodynamic calcu- 
lation offered here is asymptotically precise at small 
temperatures and, apart from missing Landau damping 
(Sec. IV Bp . just as accurate as the perturbative kinetic 
calculation in Refs. [3 5^. On the other hand, it also 
accounts for the temperature anisotropy, which is antic- 
ipated at collisionless compression or rarefaction [13] yet 
missed in Refs . liil 451. Therefore, our results complete 
those in Refs. 44 , 45j] and finally reconcile the Langmuir 



wave dynamics (particularly the temperature effects; cf. 
e.g., Refs. [H, |42|, HU) with the general action conser- 
vation theorem, which is supposed to hold for any La- 
grangian waves [To], M, H, M, M, HI . 

Besides, the above results show explicitly how the 
Langmuir wave parameters evolve. Consider, for in- 
stance, homogeneous plasma, assuming that the envelope 
shape remains fixed. Then Eq. ([6"2"]) rewrites as 

J/N = inv, (63) 

where we used that the total number of electrons is con- 
served. In the absence of Landau damping (Sec. IV Bp 
one has kvTe *C lo p , and therefore 

J w \E\ 2 /(8nuj p ). (64) 

Together with Eq. this yields (like in Ref. HH) 

E = E (N/N a ) 3 /\ (65) 

where the index denotes the initial values. Thus the 
wave field increases when the plasma is compressed and 
decreases when the plasma is rarefied. 

Finally, Eq. ([65]) also results in an effective adiabatic 
index 7 for the ponderomotive pressu re p # . To see this, 
consider the known expression for pe [69j , using that the 
field E oscillates at the frequency u> = ui p : 



PE 



\E\ 2 /(16n) 



(66) 



[In fact, Eq. ([6"B"]) itself is also derivable from Eq. ([6"3"]) , as 
shown in Appendix [Cl] Hence, from Eq. (j65p . one obtains 



PE 



\Eo\ 



Wtt \N< 



N 



V() 



3/2 



(67) 
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Therefore, for the ponderomotive pressure one has 
7 = 3/2, which is different, say, from 7 = (D + 2)/D for 
the kinetic pressure of /^-dimensional thermal electron 
gas without a wave [70l |. 



B. Wavevector. Adiabaticity conditions 



The scalings (l63)) - (l67l) hold for any wave geome- 
try, whereas the dependence of the frequency and the 
wavevector on plasma parameters may vary, as governed 
by Eq. ([4*T]) . Particularly, uj is conserved only in station- 
ary medium, and the dynamics of k is discussed below. 



1. Homogeneous plasma 

To illustrate the evolution of k, first consider plasma 
compression such that N remains homogeneous [which is 
possible at homogeneous yet not necessarily zero V • V; 
see Eq. ([54]) ]. Then the wavevector is conserved if V 
is transverse to k, an example being radial compression 
of a cylindrical plasma column with k along the axis of 
symmetry. However, if k has a component along V, the 
wavevector will evolve; specifically, 



k = ko exp 



fo<t')dt> 



with Eqs. {47]), ((49]) yielding v = v v , 

W = -k • Wk/fc 2 - V/L\ 



(68) 



(69) 



where Ly is the spatial scale on which the compression 
takes place. For instance, radial compression with V = 
%(£)r and k along V in spherical, cylindrical, and linear 
geometry equally yield vy — \. 



2. Inhomogeneous plasma 

As the next step, consider inhomogeneous plasma, for 
now assuming V = 0. In this case k can increase or 
decrease, depending on k as well as the density and 
temperature gradients, so Eq. (|68[) holds with v = ia?, 



-k-vn/k 2 



1 Pi 



(70) 



with t>ph w uip/k being the phase speed and = f2/| V£!| 
being the characteristic spatial scale. 

First, suppose that k increases. Then, on the time scale 
of order r e ~ L^/vxe, the wavelength becomes compa- 
rable to the Debye length Ad = vre/^p, regardless of 
whether the plasma inhomogeneity is due to the den- 
sity or the temperature; hence the wave decays because 
of Landau damping (see also Refs. [42|, [Zl|)- In other 
words, dissipation is negligible only at 



t < Ln/vTe 



(71) 



Thus, when compression is added, it will proceed adia- 
batically only if vyT e > 1, or 



V/L V > VTe/U 



(72) 



Assuming that the plasma average flow is entirely con- 
trolled by the large ion mass rrii 3> m e one can 
rewrite V in Eq. ([72"]) as follows. Express E from Eq. @ 
and substitute it into a similar equation for ions; hence 



9' ( V«-V-P S /M), 



(73) 



where we neglected the electron inertia and introduced 
the total kinetic pressure P E - NT [73]. Use d' t V ~ 
V 2 /L v , yielding 



V ~ c s ^Ly/Ln, 



(74) 



where c s ~ WTe\/ TO e/ m i 1S the ion sound speed. Hence 
Eq. ([72"]) rewrites as 



(75) 



Therefore, only weakly inhomogeneous plasma can be 
compressed adiabatically when k grows; otherwise a sig- 
nificant percentage of the wave energy is transformed into 
the particle thermal energy. 

Suppose now that k decreases. In this case the enve- 
lope approximation holds only on time 



(76) 



after which the wavenumber becomes zero, and thus the 
wave action is no longer conserved. Therefore adiabatic 
compression must satisfy ^ > i/yi or 

V/L v > v ph /L Q . (77) 

Assuming Eq. ([74]) . this condition hence reads as 



~ (k\n) 2 m e 



(78) 



which requires that plasma be even more homogeneous 
than in the case when k increases [cf. Eq. ([75]) ]. 



VI. CONCLUSIONS 

In this paper we show how a non-dissipative Langmuir 
wave evolves adiabatically in warm unmagnetized inho- 
mogeneous nonstationary plasma. The hydrodynamic 
calculation offered here is asymptotically precise at small 
temperatures (&Ad <SC 1) and, apart from missing Lan- 
dau damping, just as accurate as the perturbative ki- 
netic calculation in Refs. 0, |45|. On the other hand, 
it also accounts for the temperature anisotropy, which is 
anticipated at collisionless compression or rarefaction yet 
missed in Refs. J44j,j45[. Therefore, our results complete 
those in Refs. |29l. l37l l44l [4f| and finally reconcile the 
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Langmuir wave dynamics in inhomogeneous warm plas- 
mas (cf. Refs. [39|, EH ) with the general principles of 
the Lagrangian geometrical optics. 

Specifically, we derive a continuity equation [Eq. (|61[) ] 
for the wave action density, as well as the explicit condi- 
tions under which the wave action is conserved. Hence it 
is shown that, in homogeneous plasma carrying a Lang- 
muir wave, the wave field universally scales with the 
electron density as E oc A 3 / 4 . We also show that the 
wavevector evolution varies depending on the wave ge- 
ometry. Particularly, during compression k is conserved 
when aligned with the average velocity V(r,i) at homo- 



geneous density and temperature, but otherwise changes, 
with its absolute value following by Eqs. ([55|) - ([70|) . Also, 
the wave frequency to is conserved only when the plasma 
is stationary, but otherwise evolves according to Eq. (|41[) . 
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APPENDIX A: AUXILIARY VECTOR IDENTITY 

In this appendix we derive an alternative form of 

{A, B} = V • [(A • V)B + (B • V)A] - (A • V)(V • B) - (B • V)(V • A) (Al) 
for two arbitrary fields A and B. First, the expression in the square brackets above rewrites as [[zl|, Sec. 5.5-2] 

(A • V)B + (B • V)A = V(A • B) - A x (V x B) - B x (V x A); (A2) 

thus its divergence equals 

V • [(A • V)B + (B • V)A] = V 2 (A • B) - V • [A x (V x B)] - V • [B x (V x A)]. (A3) 
The chain rule for the second term on the right-hand side yields 

V • [A x (V x B)] = V • [A x (V x B)] + V • [A x (V x B )], (A4) 

where underlined are the vectors to which the differentiation by the external V applies. Because of the symmetry 
properties of the scalar triple product [[zH, Sec. 5.2-8] (here of the vectors V, A, and V x B), this also rewrites as 

V • [A x (V x B)] = (V x A) ■ (V x B) - A • V x (V x B). (A5) 

Further use that V x (V x B) = V(V • B) - V 2 B [|74|. Sec. 5.5-2]; thus Eq. (|A5|) and a symmetric expression for the 
third term on the right-hand side of Eq. (|A3|) can be put in the form 

V • [A x (V x B)] = (V x A) • (V x B) - (A • V)(V • B) + A • V 2 B, (A6) 

V • [B x (V x A)] = (V x B) • (V x A) - (B • V)(V • A) + B • V 2 A. (A7) 

Substitution of these into Eq. (|A"3|) and then Eq. (|A"3|) into Eq. (fATj) yields 

{A, B} = V 2 (A • B) - A • V 2 B - B • V 2 A - 2(V x A) • (V x B). (A8) 
In Cartesian coordinates Eq. (|A8j) finally rewrites as [jz3], Sec. 5.5-5] 

{A, B} = 2(VA j ) ■ {VBj) - 2(V x A) • (V x B), (A9) 
where summation over repeated indexes is assumed. 



APPENDIX B: PRESSURE TERMS IN THE DENSITY EQUATION 

1. General case 

In this appendix we find V • II as a function of n to the first order in e. Start off from Eq. ([9]) to get 

1 / an 3 2 P \ 

v • n - -ai + ' Wfc - a^fe J ■ ( B1 ) 
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The first term here is already of the sought form, and, since h is of order e, the second term is expressed using 

P^n(^,+P^+P^), (B2) 

as obtained from the homogeneous stationary plasma approximation [Eqs. (|T5|) , (f2T))) ]. However, for calculating the 
third term Eq. (|B2|) is not sufficiently accurate (75[, so d 2 Pjt/dxjdxi is found as follows. 
First, take d 2 /dxjdxi of Eq. ([8|), neglecting the terms of order higher than e. This yields 

d ' ( d2p ^ \ , op 5 ( V 'V) _ 



dt \dxjdxg I dxjdxe 

hnn (3^i + ^^l + ^l^l)+k J k e [ P j£ W ss + 2P js W ts + 2W js P st ) . (B3) 



The second term on the left-hand side allows an alternative representation via 

d / ~\ „ , , kjkf d 2 ( d'n 



dxjdxi V / k 2 dxjdxi \ dt 

[see Eq. ©], and the latter term in Eq. (|B4|) also equals 

d 2 ( d'n\ d' ( d' 



dxjdxi \ dt 
Hence Eq. (|B3[) can be put in the form 

d' ( d 2 P je „„ d 2 n 



{%)=U^- nlk ' k - w -' +k ' k - w - i) - m 



dt ydxjdxe ^ ^ dxjdxg J ' ( ®) 

and (J^ is given by 

5* = Znkjkt {d'tPjt) - :'./',»//(/,•, /.-Jl ., + fe^fc s W SJ ,-) + /,•,/••,(/'„ U ,, + 2P js W is + 2W js P st ). (B8) 
Using the slow component of Eq. ([3|) in the form 

d' t P je = -PjiWss - P 3S W ls - P sl W 38 (B9) 

and also Eq. (|B2[) , rewrite 8^ as 

= -4nk • (WGPk), G = i - kk//c 2 . (BIO) 
Since $ is a rapidly oscillating function with a slow envelope of order e, Eq. (|B6[) is integrated as 

a 2 P,£ _ 3p <9 2 n + 

dxjdxi ^ dxjdxi 

Substitute Eq. (IBlljl for the third term in Eq. (|Blj) . together with Eq. (jB2|) for the second term; then one gets 

v ■ n =~ c »i£k - *- («<■ + + f (I kW - h ) • 66k (B12) 

Since the wave is assumed propagating along a local principal axis of the temperature tensor [64j |. one has GCk = 
Gk x const. Yet Gk = 0, so Eq. (|B12[) is simplihed, and, using inki m dn/dxi, one finally obtains 



kike \ dCef d 



2. Isotropic temperature 



For the isotropic temperature case, Eq. (|B13[) gives [62| 
V 11 = -34 e V 2 n-6Vn- Vu| e , (B14) 



which as well can be obtained by substituting Eq. (80) 
in Sec. IV.3 of Ref. [I| into our Eq. (EE]). 
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Alternatively, Eq. (|B14|) can be derived using a phe - 
nomenological adiabatic law (cf., e.g., Ref. [7q , Ref. []z3, 
Sec. 5.1], Ref. [[E3], Sec. 3.5]) 

(d t +V e -V)( Pe N-^)=Q, (B15) 

with p e being the electron scalar pressure, including the 
slow and the quiver parts: p e = p + p- Here 7 = 3 (corre- 
sponding to one-dimensional adiabatic oscillations [7(3]) 
is an extrapolation from the homogeneous plasma case, 
for which the exact solution is known from a more rig- 
orous Iwdrodynamic treatment, like in our Sec. IIIII or 
Refs. }48l l49l. pQ ] . or the complete kinetic treatment [[52|. 
Chap. 8], [531 ]. To show this, introduce the plasma ele- 
ment Lagrangian displacement £ such that [u], [78|, [7§] 

N + V • (£N) = 0, (B16) 
p + £ • Vp = m e a 2 (N + £ ■ VAT), (B17) 

where a 2 = 3u| e . Then, from Eq. © with P = pi and 
P = pi, one gets 

n = -V(na 2 ) 



[nq+V(£-q)]/iV. 



(B18) 



Here q = a 2 VA — Vp/m e is of order e, and therefore one 
can take d' t £ w V (cf., e.g., Eq. (4 .6) in Ref. [z|), so 
£ ps mk/A; 2 , as flows from Eq. (|B16[) . Hence 



n w -V(na 2 ) - nAf-^q, 

yielding 

V n « -V 2 (na 2 ) - ircJV^kG • q. 
Using that kG = 0, one finally obtains 



V n w -a 2 V 2 n - 2Vn • Va 2 . 



which is equivalent to Eq. (|B14|) . 



(B19) 



(B20) 



(B21) 



APPENDIX C: PONDEROMOTIVE PRESSURE 

The conservation of the Langmuir wave action allows 
to calculate the effective stress tensor due to the wave, 
which is done as follows (see also Refs. [UE^, 80]). For 
simplicity suppose homogeneous cold stationary plasma 
volume V and assume that it is adiabatically deformed 
as defined by an infinitesimal displacement field £(r), re- 
sulting in the strain tensor 



w = [(V£) + (V£) T ]/2. 



Hence the stress tensor <r is found [8l| 



(CI) 



i as 

V owj£ 

where (£ = VJuj is the wave total energy inside V. Be- 
cause V J is conserved, Eq. ()C2|) rewrites as 



Juj p dN 



2N dw 



(C3) 



where we used that u> — u> p (N). The density perturba- 
tion due to the strain is SN = -NV ■ £ [cf. Eq. (|BT6|) ]. 
Since V • ^ = w ss , and dw ss /dwjg — 5jj>, this yields 



-(Jwp/2) i. 



(C4) 



Therefore the stress due to the wave field is isotropic and 
appears as an effective pressure pe = Joj p /2. (Thermal 
correction would also yield an anisotropic component to 
the wave stress tensor [H3,[6§].) Using Eq. (l64|) . one then 
recovers the expression [Eq. ([66]) ] for the ponderomotive 
pressure in cold plasma carrying a field E which oscillates 
at the frequency lu — ui p [69( . 
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